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We develop and analyze an ordinary differential equation model to investigate the transmission dynamics of 
releasing Wolbachia-miecteA mosquitoes to establish an endemic infection in a population of wild uninfected 
mosquitoes. Wolbachia is a genus of endosymbiotic bacteria that can infect mosquitoes and reduce their ability to 
transmit dengue virus. Although the bacterium is transmitted vertically from infected mothers to their offspring, 
it can be difficult to establish an endemic infection in a wild mosquito population. Our transmission model for 
the adult and aquatic-stage mosquitoes takes into account Wolbachia-indnced fitness change and cytoplasmic 
incompatibility. We show that, for a wide range of realistic parameter values, the basic reproduction number, Rq, 
is less than one. Hence, the epidemic will die out if only a few Wolbachia-infected mosquitoes are introduced into 
the wild population. Even though the basic reproduction number is less than one, an endemic Wolbachia infection 
can be established if a sufficient number of infected mosquitoes are released. This threshold effect is created by a 
backward bifurcation with three coexisting equilibria: a stable zero-infection equilibrium, an intermediate-infection 
unstable endemic equilibrium, and a high-infection stable endemic equilibrium. We analyze the impact of reducing 
the wild mosquito population before introducing the infected mosquitoes and observed that the most effective 
approach to establish the infection in the wild is based on reducing mosquitoes in both the adult and aquatic 
stages. 

Keywords: dengue; vertical transmission; backward bifurcation; cytoplasmic incompatibility; basic reproduction 
number 


1. Introduction 

We use a disease transmission model to investigate the conditions for releasing Wolbachia-iniected 
mosquitoes that will establish an endemic infection in a population of wild uninfected mosquitoes. The 
Wolbachia infected mosquitoes are less able to transmit dengue virus, and the goal of the modeling effort 
is to better understand how this bacterium can be used to control vector-borne diseases. We observed that 
the basic reproductive number for the Wolbachia model is less than one for typical model parameters, and 
so small Wolbachia infestations will die out. However, the model predicts that there is a critical threshold, 
and if a sufficient number of infected mosquitoes are released, then an endemic Wolbachia-infected popu¬ 
lation of mosquitoes can be established. We also investigated how this critical threshold can be reduced by 
first decreasing the population of wild mosquitoes. 

Dengue is the world’s most significant and widespread arthropod-borne viral disease [l|. Each year, 400 
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million people are infected with dengue virus in more than 100 countries while the other one-third of 
the world’s population is at risk. To date, there are no vaccines or specific therapy available. Traditional 
control strategies focusing on reducing population of Aedes mosquito vectors have failed to slow the current 
dengue pandemic, especially in tropical communities Q]. The main vector, Aedes aegypti, rebounds in many 
areas and the secondary vector, Aedes albopictus keeps expanding its geographic distribution, leading to 
30 fold increase in cases over the past 50 years 0], which necessitates effective novel alternatives to break 
dengue transmission cycles targeting Aedes aegypti and Aedes albopictus. 

Increasing attention has been paid to controlling the spread of dengue by targeting mosquito longevity by 
introducing genetically modified mosquitoes or introducing endosymbiotic Wolbachia bacteria to shorten 
the mosquito lifespan [3, S, [^. That is, Wolbachia-vaiecied mosquitoes are released to create a sustained 
infection in the wild (uninfected) population. If the infection is sustained, then the wild infected mosquitoes 
will be less effective in transmitting dengue fever. We create and analyze a mathematical model to help 
understand the underlying dynamics of Wolbachia-iTdected mosquitoes that are needed to create a sustained 
endemic Wolbachia infection. Once the Wolbachia disease transmission model is well understood, one of 
our future goals will be to couple this model with a mosquito-human model for the spread of dengue. 


1 . 1 . Wolbachia Bacteria 


The wMel strain of Wolbachia pipientis bacteria is a maternally inherited endosymbiont infecting more 
than 60% of all insect species. This strain has the ability to alternate host reproduction through partheno¬ 
genesis, which results in the development of unfertilized eggs, male killing, feminization, and cytoplasmic 
incompatibility (Cl) [3, that prevents the eggs from forming viable offspring. The latter includes strate¬ 
gies for both the suppression and replacement of medically important mosquito populations. Cytoplasmic 
incompatibility is an incompatibility between the sperms and eggs induced by Wolbachia infection and 
has received considerable attention as a method to control vector-borne diseases (lOtl. Wolbachia even in¬ 


duces resistance to dengue virus in Aedes aegypti 11[ and limits transmission of dengue virus in Aedes 
albopictus (l^ . 

Uninfected females only mate successfully with uninfected males, while infected females can mate success¬ 
fully with both uninfected and infected males 1^. If a male fertilizes a female harboring the same type of 
infection, the offspring still can survive [^. When Wolbachia-infected males mate with uninfected females, 
or females infected with a different Wolbachia strain, then the Cl often results in killing the embryos. Q|. 
Therefore, infected females have a reproduction advantage over uninfected females due to protection from 
Cl Q. 

To successfully transmit dengue virus, a vector must imbibe virus particles during blood-feeding and 
survive to the point that the pathogen can be biologically transmitted to the next vertebrate host [l^. This 
time period, called extrinsic incubation period (EIP), varies with ambient temperature, many other climatic 


factors, and characteristics of the vector-parasite system 1^. Vectors that survive long enough to transmit 




the pathogen are called effective vectors Il5|. Typically dengue virus has an incubation period as long as 


two weeks to transmit through Aedes aegypti populations [1[. A life-shortening strain of Wolbachia may 




halve the life span of Aedes aegypti |l7l |. Wolbachia infection may reduce the rate of disease transmission 
due to the reduction on the lifespan of infected mosquitoes or the interference with mosquito susceptibility 
to dengue virus. 
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1 . 2 . Existing Mosquito-Wolbachia Models 


Ordinary differential equation (ODE) models have been developed to explore key factors that determine 
the success of applying Wolbachia to dengue control. A single-sex model for Wolbachia infection with both 
age-structured and unstructured models were presented to study the stability and equilibrium based on 
the assumption that Wolbachia infection leads to increased mortality or reduced birth rate 1^. A model 
assuming a fixed ratio of females and males addressed how pathogen protection affects Wolbachia invasion 


19|. Age-structured and unstructured models combining males and females were found to be different in 


terms of existence and stability of equilibrium solutions . A stochastic model for female mosquitoes was 
developed to investigate the impact of introduction frequency on establishment of Wolbachia 20|. 

Discrete-time models explored the impact of the type of immigration and the temporal dynamics of 
the host population on the spread of Wolbachia^ assuming equal sex ratio between males and females 


21[. Discrete generation models for female mosquitoes were built to understand unstable equilibrium 


produced by reduced lifespan or lengthened development 2^. Reaction-diffusion and integro-difference 
equation model has been used to analyze the impact of insect dispersal and infection spread on invasion 
of Wolbachia 


23|. 


An ordinary differential equation (ODE) model was developed to evaluate the desirable properties of 
the Wolbachia strain to be introduced to female mosquitoes, assuming that Wolbachia-infected mosquitoes 
have reduced lifespan and reduced capability to transmit dengue, and equal fraction of male and female 
mosquitoes 2^. A continuous time non-spatial model and a reaction-diffusion model incorporating lifespan 
shortening and Cl were developed to study factors that determine the spatial spread of Wolbachia through 
a population of female Aedes aegypti mosquitoes assuming constant population size and perfect maternal 
transmission of Wolbachia 2^. A two-sex deterministic model with deterministic immature life stages and 


stochastic female adult life stage was developed to understand Wolbachia invasion into uninfected host 
population (26|. 

A single strain model, two strain model, and spatial model were developed to study whether multi-stain 
of Wolbachia can coexist in a spatial context 27|]. A two-sex ODE model taking into account different 
death rates, but the same egg laying rates of Wolbachia-iniecied and Wolbachia uninfected mosquitoes 


28l | showed the basic reproduction number is always less than one, and the complete infection equilibrium 


is locally asymptotically stable (LAS) due to positive determinant of the Jacobian matrix for the system. 
Simulations showed that dengue epidemics will not occur when Wolbachia infection is sufficiently prevalent 


28l |. A two-sex mosquito model assuming complete vertical transmission and equal death rates for male 


and female mosquitoes was developed and four steady states were found [29|. 

Most of these models consider either a single-sex model for adult mosquitoes, or assume a fixed ratio 
between the number of male and female mosquitoes. Also, most of the models assume homogeneous death 
rates and egg laying rates for Wolbachia-free and Wolbachia-mfected mosquitoes. Our model addresses 
both of these issues. Our emphasis is to understand how Wolbachia infection can be established in a wild 
population of mosquitoes. 


1 . 3 . Results 

We proposed a compartmental two-sex model to investigate the underlying mechanisms that may contribute 
to invasion and sustainable establishment of Wolbachia in mosquito populations. We assigned female and 
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male mosquitoes to different classes to understand corresponding roles that they are playing in the spread 
of Wolbachia in mosquito populations. 

We showed that, for a wide range of realistic parameter values, the basic reproduction number, Ro, for 
this model is less than one. Hence, the epidemic will die out if only a few Wolbachia-infected mosquitoes 
are introduced into the wild population. Even though the basic reproduction number is less than one, an 
endemic Wolbachia infection can be established if a sufficient number of infected mosquitoes are released. 
This threshold effect can be explained as a backward bifurcation with three coexisting equilibria: a stable 
zero-infection equilibrium, an intermediate-infection unstable endemic equilibrium, and a high-infection 
stable endemic equilibrium. 

If the number of infected individuals is below the unstable endemic equilibrium, then the infection decays 
to the zero-infection equilibrium. Conversely, if the number of infected mosquitoes is greater than unstable 
endemic equilibrium, then the solution tends to the stable high-infection equilibrium. We identified the rela¬ 
tionships between dimensionless combinations of model parameters and the initial conditions for Wolbachia 
to be attracted to the high-infection state. As expected, the number of infected female mosquitoes needed 
to be released to establish the infection in a wild Wolbachia-iree population decreases as Rq increases to 
one. We analyzed the impact of reducing the wild mosquito population before introducing the infected 
mosquitoes. We found that the most effective approach of reducing the number of infected mosquitoes 
needed to establish a wild Wolbachia-iniected population requires reducing wild mosquito populations in 
both the adult and aquatic stages before the release. This could be accomplished by recursive spraying, or 
a combination of spraying and larvae control. 

Our main findings are: 

(1) Three equilibria, disease free equilibrium (DEE), endemic equilibrium (EE), and complete infection 
equilibrium (CIE) coexist when < 1- Disease free equilibrium is a steady state when all individuals 
are Wolbachia-iree. Endemic equilibrium is a steady state when some individuals are Wolbachia-iree, 
the rest are infected with Wolbachia. Complete infection equilibrium is a steady state when all indi¬ 
viduals are infected with Wolbachia. 

(2) The backward bifurcation analysis of our Wolbachia transmission model predicts that if Rq < 1, then 
there is a critical threshold for the number of infected mosquitoes released in the wild before the 
infection can be established. If we release too few Wolbachia infected mosquitoes, then Wolbachia 
infection will die out. 

(3) Killing both aquatic state (eggs and larvae) and adult mosquitoes before releasing the infected 
mosquitoes greatly increases the chance that the infection will be established. Our model quantifies 
what fraction of wild mosquitoes must be killed before releasing Wolbachia infected mosquitoes. 

After introducing the mathematical model, we summarize the key results from the analysis and numerical 
simulations. We conclude with a discussion of the relevance, importance, and future directions for this work. 


2. Description of Model Framework 

We developed an ODE model incorporating adult females (E), adult males (M), and an aggregated aquatic 
(A) stage that includes the egg, larvae, and pupae stages. The population dynamics of mosquitoes without 
taking into account Wolbachia is in the Appendix, Equation [9l The vertical transmission of Wolbachia 
from infected females to their offspring is a key factor in establishing an endemic infected population. 
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Mosquitoes are grouped into six compartments: susceptible aquatic stage, infected aquatic stage, ^4^, 
susceptible female mosquitoes, Fu, infected female mosquitoes, susceptible male mosquitoes, M„, and 
infected male mosquitoes, M^. The eclosion rates of susceptible female and male mosquitoes hatching from 
eggs are and 'il)(l—6)Au, respectively. Similarly, the birth rates of infected female and male mosquitoes 

are TpOA^ and ?/>(l — 6)A^. Death rates of uninfected male mosquitoes and infected male mosquitoes are 
IXaAu and fJ.aAw Death rates of uninfected female mosquitoes and infected female mosquitoes are HjuFu 
and ^fyjFw Death rates of uninfected male mosquitoes and infected male mosquitoes are HmuMu and 
HmwAIw Development rates of uninfected aquatic stage and infected aquatic stage of mosquitoes are 
and ijjAu]. The model parameters are described in Tabled) 

The model (Fig. d]) describing population dynamics of aquatic stage, adult male, and adult female 
mosquitoes is given by: 


dAu 

dt 


dt 

dFu 

dt 

dFy, 

dt 

dMu 

dt 

dMyj 

dt 


Fuu T WU + B 

WW ) - fJ,aAu - i^Au 

V^{B 

WU + B 

WW ) f^a^w '^^w 

bflpAy^ fJ>fuFu 

bf'lpA^ 

bm 

bm 'IpAru^ 


(la) 

(lb) 

(l c) 

(l d) 

(le) 

(l f) 


Because the vertical transmission and birth rates, B**, depend on the sex of the infected or uninfected 
parents, the model included the four egg laying situations 


Btini — (biiF^iTTL^, { 1 


Bijnn - 0 


K„ 


Bwu - 4^wFu!l^U I 1 


Ka 

Bww — fpwFwTTlw ( 1 „ 


(2a) 

(2b) 

(2c) 

(2d) 


Here and rriyj = 1 — are the fractions of uninfected and infected male mosquitoes. B^u 

is the egg laying rate of uninfected females mating with uninfected males. B^w is the egg laying rate of 
uninfected females mating with infected males. B^u is the egg laying rate of infected females mating with 
uninfected males. Na is the total number of aquatic stage of mosquitoes, and Ka is carrying capacity of 
aquatic stage of mosquitoes. These equations reflect the observations in Table [2] that 

• Mating of uninfected males with uninfected females produce uninfected offspring. 

• Mating of infected males with uninfected females leads to death of embryos before hatching due to 
cytoplasmic incompatibility. 

• Uninfected males and infected females produce a fraction, denoted by Vw, of infected offspring by vertical 
transmission. 
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Figure 1. The birthing rates m capture that when the uninfected males mate with uninfected females, they produce uninfected offspring. 
When infected males mate with uninfected females, then Cl causes the embryos to die before hatching. Uninfected males mating with 
infected females produce a fraction, denoted by Vw^ of infected offspring by vertical transmission. Cross of infected males with infected 
females produces a fraction of infected offspring. 


• Cross of infected males with infected females produces a fraction of infected offspring. 


2.1. Model Analysis 

We compute the basic reproduction number and the equilibria, and analyze the stability of the equilibrium 
points. 

We use the next generation matrix approach to compute the basic reproduction number 30|]. Only 
infected compartments are considered for ease of computation: 


dt 


1 T 


Fu, Mu 


= ^-y = 


{VwFu,(pw^u “b '^wFu,4^w'^w^ 

(i-fe)' 


l^a-^w “1“ 

bfipAu, 


— 

f^fw^w 

hmll^Au, 



l^mw^w 


where ^ = {^i) is a vector for new infected, and F = {Fi) is a vector for transfer between compartments. 
Jacobian matrices for transmission, F, and transition, V, 30| are defined as: 


d^i{x^) 

, V = 

\dyi{x^)] 

[ dxj \ 


[ dxj \ 


(4) 


where represents the disease free equilibrium, and Xj is the number or proportion of infected individuals 
in compartment j, j = 1,2, • • • , m. 
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Mu'. Number of susceptible male mosquitoes 

M^: Number of male mosquitoes infected with Wolbachia 

Fu'. Number of susceptible female mosquitoes 

Fw'. Number of female mosquitoes infected with Wolbachia 

Au'. Number of susceptible aquatic stage of mosquitoes 

A^: Number of infected aquatic stage of mosquitoes 

Na- Total number of aquatic stage of mosquitoes 

Ka‘. Carrying capacity of aquatic stage of mosquitoes 

bj: Fraction of births that are female mosquitoes. 

b^: Fraction of births that are male mosquitoes = 1 — bf. 

rriw- Fraction of the male mosquitoes that are infected = M^/(Mu, + M„). 

rriu- Fraction of the male mosquitoes that are uninfected = 1 — m^. 

Vw- Fraction of infected mosquito eggs produced by infected female mosquitoes. 

Vu- Fraction of uninfected mosquito eggs produced by infected female mosquitoes = 1 — Vw 
(pu- Per capita egg laying rate by Wolbachia-hee mosquito eggs. Number of eggs/time 

Per capita egg laying rate by Wolbachia-iniecied mosquito eggs. Number of eggs/time 
'i/j: Per capita development rate of mosquito eggs. Time“^ 

fia- Per capita death rate of aquatic stage of mosquitoes. Time“^ 

fifu- Per capita death rate of uninfected female mosquitoes. Time“^ 
fif^: Per capita death rate of infected female mosquitoes. Time“^ 

lamu- Per capita death rate of uninfected male mosquitoes. Time“^ 

IJ-mw- Per capita death rate of infected male mosquitoes. Time“^ 

Table 1. State variables and parameters for the model Q 


Male 

Female 

Offspring 

Uninfected 

Uninfected 

Uninfected 

Uninfected 

Infected 

Some offspring are infected 

Infected 

Uninfected 

Unviable 

Infected 

Infected 

Some offspring are infected 


Table 2. The model allows for Wolbachia to be transmitted vertically from infected parents to their offspring. The offspring of male and 
female uninfected mosquitoes are infected. Some of the offspring of male and female infected mosquitoes are infected, as are the offspring 
of an uninfected male and an infected female mosquitoes. The offspring of an infected male and uninfected female mosquito are unviable. 


The unique disease free equilibrium is 


A'i = Ka{l- 


Rqu 


F^u = bi 




A?, 


hfu 

rO t. ^ 


M" = 6, 


l^mu 


■A^,. 


Ai = F^ = Mi:, = 0 , 


?o 
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where Rqu = is the threshold for lf^oZ6ac/iM-free offspring, is the total number of eggs 

laid by uninfected female mosquitoes, -^-r is the probability that aquatic stage of mosquitoes survive to 

^a~rW 

the point when they develop into adult mosquitoes. Rqu is the number of female eggs that develop into 
adult mosquitoes. When > Ij then the mosquito population may grow; otherwise, the population will 
decrease. 

The Jacobian matrix of ^ evaluated at DFE is: 


0 

Vwmlct)^ (^1 - 

o' 


0 

j;»(Ala+b)M/^b™ 

o' 



bfif} 

0 

0 


bfi) 

0 

0 

bm'tp 

0 

0. 


bm'ip 

0 

0_ 


where ml = m^+m° • 

The Jacobian matrix of "Y evaluated at DFE is: 


V = 


Ma + t/’ 
0 
0 


0 

0 


0 

0 


The next generation matrix for infected compartments at disease free equilibrium is: 


FV-^ 


0 

Ma+b 

bm'ip 

-TT+b 


0 

0 


0 

0 

0 


The basic reproduction number is the largest abosulute eigenvalue of FV denoted by p{FV ^): 


Ro = p{FV-^) 



(5) 


Wolhachia-iniecied and Wolhachia-ixQQ and females produce and eggs during their lifetime, 

respectively. Hence, (/^) ratio of the number of eggs produced by Wolbachia-mfected females 

to the number of eggs produced by lTo/6ac/iia-free females during their lifetime. Rq is geometric mean of 
vertical transmission rate and the ratio of the total number of eggs produced by Wolbachia-iniected females 
to the total number of eggs produced by Wo/tec/iia-free females. 

When the variables are ordered as: [A^ F^ F^j M^]^, Jacobian matrix at equilibrium 

[H* FI M* AIj for system of equations ([I]) is: 


c 

a 

b 

C2 

02 

h2 1 




bftp 

~Pfu 

0 

0 

0 

0 




bm'4^ 

0 

l^mu 

0 

0 

0 


'a 

C 

C3 

0 

h 

Cl 

Ol 

hi 


D 

B 

0 

0 

0 

bfijj 

l^fw 

0 




0 

0 

0 

bmi’ 

0 

l^mw_ 
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where 


A* A* 

a = ( 1 - “ ^ 


Kn 


b = (t)uml 


M* + M* 


1 - 


A* 4- A* 
Ka 


C = -(/ia + V’) - - (pwVu^ 

_ K + K 

I 1 

bi = 0 

T7 

Cl = -(//a + "0) - 4>wVn 


K, 


CL2 — 4‘w'^u ( 1 




Ka 


b2 = -(t>uml 


p* 


M* + M* 


A* _1_ /I* 
^ _ -^u ' -^w 


Ka 


P* p* 

1 ■¥ U 1 W 

C2 - 4^u'^u 4^w^u 

b3 = 0 


C3 — (pw'^'U 


Ka 


At disease free equilibrium, D = 0, then the eigenvalues of J are eigenvalues of A and B. Matrices A and 
B are: 



(Mo + V’) ~ 

ll/u(lla+b) 

0 


-{pa + '4’) 


0 


bfi!) 



A = 

bf-lP 

~Tfu 

0 

, B = 

bflP 

^fw 

0 


bmi’ 

0 

f^mu_ 


bm'4^ 

0 

t^mw_ 


If Rou > 1, then all eigenvalues of A are negative. If i?o < Ij then all eigenvalues of B are negative. 
Therefore, the system ([T|) at disease free equilibrium is LAS whenever Rqu > 1 and i?o < 1- 


2.1.1. Complete Vertical Transmission 

If i?o < 1 and Rqw > 1 and the vertical transmission is 100%, {vw = 1), then the ratio of the infected to 
the uninfected aquatic stage mosquitoes, k, is 


A 


t^mw / 72)—2 


k = ^ = ^{Ro^-l) , 

■^U f^TTlU 


and there is a unique endemic equilibrium: 


“ l + k 

4* _ L4* 
■^w 


(1-^oi) 


(6a) 

(6b) 
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/?* 

u 

/?* 

w 

M* 

M^. 


h ^ A* 
fJ-fu 

kF* 


if; 


■A* 


f^mu 

kM:, 


(6c) 

(6d) 

(6e) 

(6f) 


where Rn — 

wneie ito^j — • 

If we further assume that Hmu 

(CIE) is: 


fimw, then k = R 


-2 

0 


1 and the unique complete infection equilibrium 


Au = 0 

Aw = Ka{l- i?oJ) 

Fu = 0 

Fw — bf A-w 
f^fw 

Mu = 0 

xb 

AIw — bra Aw 
f^mw 

Jacobian matrix of system of equations (H]) at complete Wolbachia infection equilibrium for = 1 is: 


-Ra - Ip 

0 

0 1 

0 

0 

0 




bfXP 

~Rfu 

0 1 

0 

0 

0 




bmp’ 

0 

l^mu 1 

0 

0 

0 


'a 

0 

bfip{l-R'^)<t)^ /, p-l\ 

0 

0 1 


(/Xa+V’lM/™ 

0 


c 

B 

0 

0 

0 1 

bfip 

l^fw 

0 




0 

0 

0 1 

brapJ 

0 

l^mw 





Eigenvalues of Jew are composed of three eigenvalues of A and three eigenvalues of B. The eigenvalues of 
A are all negative. Characteristic polynomial of B is: 


\2 . I I 4'wbf'll^ 

A + ( Rfw + — - 

l^fw 


^ A + 4>wbfi> - ina + 'llj)bl‘fw = 0 


(7) 


If Rqw > 1, all eigenvalues of B are negative. Therefore, the complete infection equilibrium is LAS whenever 


Row > 1, as shown in Figure 2(a) 
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Parameter 

Baseline Value 

Ranges 

References 

bf 

0.5 

0.34 - 0.6 

\2l-M 

4^u 

50/day 

0-75 

[32] 

4^w 

51/day 

0-75 

[32] 

tp 

0.01/day 


Assume 

Aia 

0.02/day 


Assume 


0.9 

0-1 

Assume 

k-fu 

0.061/day 

1/55 - 1/11 

[32, 35] 

l^fw 

0.068/day 

1/55 - 1/11 

[32, 35] 

l^mu 

0.068/day 

1/31 - /7 

[32,M, 35] 

l^mw 

0.068/day 

1/31 - /7 

[^, 3^ 

Ka 

100,000 


Assume 


Table 3. Baseline values for parameters in Model Q. 


Vertical 

Transmission 

Disease Free 
Equilibrium (DEE) 

Endemic 

Equilibrium (EE) 

Complete Infection 
Equilibrium (CIE) 

Complete 

Ro < 1 and Rqu > 1, 
LAS 

Ro < 1 and Rqw > Ij unstable 

Rou, > 1, LAS 

Incomplete 

Ro < 1 and Rqu > 1, 
LAS 

0.5 < < 1 and Avwivw — 1) < 

Ro < 1, when k = ki, LAS, when 
k = k 2 , unstable 

does not exist 


Table 4. Threshold condition for existence of disease free equilibrium, endemic equilibrium, and complete infection equilibrium and 
their stability. Rqu is the threshold for Wolbachia-iree mosquito population, and Rqw is the threshold for Wolbachia-mfected mosquito 
population. Only when Rqu > 1, Wolbachia-iree may grow, and only when Row > 1, Wolbachia-iniected population may grow. 


2.1.2. Incomplete Vertical Transmission 

When vertical transmission is incomplete, i.e., 0 < Vw < 1, then at endemic equilibrium, 



( 8 ) 


We 


assume Hmu = IJ-mw, and let A: = ^. When (vu^ - t’lo) 


-2 

0 


2y/VwVu. Equation [7] has two roots: ki = 

If 0.5 < < 1 and ‘ivwivw — 1) < Rq < 1, two endemic equilibria exist as shown in Figure 2(b), 2(c) 


> 4vuVw{Ro‘^ “ 
and k 2 = 


that is, i?o > 


and |2(d)1 When A: = A:i, the endemic equilibrium is LAS. When k = k2, the equilibrium is not LAS and 
backward bifurcation occurs. If < 0.5 and Rq < 1, then endemic equilibrium does not exist, only DEE 
exists. When ki = k 2 , then Rq = 2y/VwVu- Let Rq = 2y/VwVu, which is the intersection of unstable and 
stable endemic equilibrium. When Rq > 1, a unique endemic equilibrium point exists with A: = fei. It is 
proven to be LAS by numerical simulations. 
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(d) v^u = 0.75. 


Figure 2. Bifurcation diagrams for Wolbachia vertical transmission, (pu and are varying, other parameter values are the same as 
those baseline values in Table 0 Denote the intersection of two endemic equilibrium, that is, the intersection of the black dashed line 
and x-axis, as Rq = y/AvwVu- When Rq < 1 and Vw < 0.5, no endemic equilibria exist. When Rq < 1 and Vw > 0.5, as the vertical 
transmission rate increases so does Rq and the LAS equilibrium approaches a constant. If we increase the number of infected females, then 
the endemic equilibrium may become stable endemic or complete infection equilibrium. If we decrease the number of infected females 
at endemic equilibrium, then the endemic equilibrium may become disease free equilibrium. WIF denotes Wolbachia-infected female 
mosquitoes. 


3. Results 


We observed that Wolbachia can persist when Rq < Rq, where Rq is the turning point of the backward 
bifurcation. Three equilibria, namely, disease free equilibrium, endemic equilibrium, and complete infection 


equilibrium coexist when Rq < 1, Rq^ > 1, = 1, and Rqu > 1 as shown in Figure 2(a) The disease free 


equilibrium is LAS whenever Rq < 1 and Rqu > 1; and complete infection equilibrium exists and is LAS as 
long as Rqw > 1- The unique endemic equilibrium is not LAS, which can become disease free equilibrium by 
decreasing the number of infected individuals or become complete infection equilibrium by increasing the 
number of infected individuals. Figures 2(b) 2(c), and |2(d)] showed that two endemic equilibrium points 
exist when < Rq < 1 and Vw > 0.5, and only one of them is LAS proven by numerical simulations. 

When Vw is larger, R* is closer to one. The condition for the existence of the equilibria and their stability 
are summarized in Table 01 
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Figure 3. Thresholds for fraction of infected individuals vary with reproduction number. AuO+Aujq = A®, FuO + FwO = F®, MuO+M-wO = 
M^. When Rq < 1, the smaller Rq is, the larger number of infected female mosquitoes are needed to be released for Wolbachia to be 
endemic. The Wolbachia infection is only sustained if the fraction of WIF mosquitoes is above the red dotted line. 


Initial condition thresholds for an epidemic to occur vary with the vertical transmission rate as shown 
in Table [5l When Vw = 0.9, the epidemic will spread if at least 40% of the population are initially infected. 
When Vw = 0.95, the epidemic will spread if at least 34% of the population are initially infected. When 
= 1, the epidemic will spread if at least 28% of the population are initially infected. When the vertical 
transmission rate is high, the threshold for the number of initially infected individuals to start Wolbachia 
epidemic is low. 

The initial condition threshold for an epidemic to occur varies with the ratio of death rates of Wol- 
bachia-miected male mosquitoes to death rates of Wolhachia-ixee male mosquitoes, while hxing other 
parameters. The larger is, the larger the initial number of infected individuals required to start a Wol- 
bachia epidemic. With higher vertical transmission rate, a smaller percentage of Wolbachia carriers can 
invade, as shown in Table El Similarly, if we increase (pw or fifu, or decrease (pu or /r/iu, then the threshold 
for initial infection is smaller. Increasing ^ will increase the epidemic threshold for initial infection. 
The thresholds for fraction of initial infections decreasing with the basic reproduction number is shown in 

Figure El 

The reproduction number is very sensitive to the vertical transmission rate, egg laying rates of Wolbachia- 
infected mosquitoes, egg laying rates of Wolbachia-fiee mosquitoes, and death rates of Wolbachia-iniected 
female mosquitoes and Wolbachia-hee female mosquitoes as shown in Equation ([5]). The reproduction 
number varies directly with either the vertical transmission rate, the egg laying rates of Wolbachia-infected 
mosquitoes, or the death rates of Wolbachia-free female mosquitoes. The reproduction number varies 
inversely with egg laying rates of Wolbachia-free mosquitoes, or the death rates of Wolbachia-infected 
female mosquitoes, but towards opposite direction. 

We compared five strategies before the release of Wolbachia-iniected female mosquitoes: 

• DFE: releasing Wolbachia-inlecied female mosquitoes at the disease free equilibrium, 

• KHA: first killing half of the aquatic stage of mosquitoes, 

• KHM: first killing half of the wild adult mosquitoes, 

• KHM2: first killing half of the wild adult mosquitoes, and then killing half of the adult mosquitoes again 
after two weeks, 

• KHMA: first killing half of the wild mosquitoes and half aquatic stage of mosquitoes. 
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Initial condition 

Vw = 0.9 

Vw = 0.95 

Vm — 1 

^uO 

< 60%A° 

< m%Al 

< 72%A0 

■^wO 

> 40%AO 

> MVoAl 

> 28%A0 

Fuo 

< 60%F0 

< 66%F^ 

< 72%F^ 

^wO 

> 40%F0 

> 34%F0 

> 28%F^ 

Muo 

< 60%M0 

< 66%M0 

< 72%M0 

AIwo 

> 40%M0 

> 34%M0 

> 28%M0 


Table 5. Initial condition thresholds for epidemic to occur with different vertical transmission rates. When the number is Wolbachia- 
carrying mosquitoes is above the threshold, Wolbachia-'miected mosquitoes can petsit, otherwise, they are wiped out by Wolbachia-iree 
mosquitoes. 


Scenario 

minimum release ra¬ 
tio when Rq = 0.85 

minimum release ra¬ 
tio when Rq = 0.9 

minimum release ra¬ 
tio when Rq = 0.95 

DFE 

1.007 

0.571 

0.259 

KHA 

0.549 

0.317 

0.151 

KHM 

0.478 

0.244 

0.107 

KHM2 

0.329 

0.165 

0.072 

KHMA 

0.210 

0.112 

0.051 


Table 6. Different population suppression strategies applied before release of Wolbachia infected female mosquitoes can reduce the 
minimum number of Wolbachia-lnfected mosquitoes that can lead to persistence of Wolbachia. The top row indicates that you have to 
release 1.007 times as many infected female mosquitoes as there are wild infect mosquitoes to establish an infection in the wild that starts 
at the DFE when Rq = 0.85. The number of infected mosquitoes that need to be released to establish an infection decreases with the 
strategies: KHA denotes killing half of the aquatic stage mosquitoes, KHM denotes killing half of adult wild mosquitoes, KHM2 denotes 
killing half adult of wild mosquitoes once, then kill half adult wild mosquitoes again after two weeks, and finally KHMA denotes killing 
half adult wild mosquitoes and half aquatic stage of mosquitoes. 


The ratios of the minimum number of Wolbachia-infected female mosquitoes that can lead to persistence of 
Wolbachia to the number of female mosquitoes at disease free equilibrium are listed in Table [6] in decreasing 
order. Notice that killing the adult mosquitoes only once is not an effective strategy because the aquatic 
stage mosquitoes hatch and quickly replace the wild uninfected population. Killing both the adult and 
aquatic (larvae) stage mosquitoes before releasing the infected mosquitoes is the most effective strategy. 


4. Discussion 

We developed a model considering two sex of aquatic stage and adult mosquitoes, diversity in the death 
rates of Wolhachia-miecied mosquitoes and Wolhachia-iiee mosquitoes, and egg laying rates of Wolbachia- 
'miected female mosquitoes and Wolbachia-fiee female mosquitoes. The general model is not constrained to 
particular weather condition, specific Wolbachia strains, or specific mosquito species, and it can be easily 
adapted to Aedes aegypti or Aedes albopictus at any location with parameters calibrated using realistic 
environmental factors, such as temperature and rainfall etc. 

We found conditions for the existence of multiple equilibria and backward bifurcation. If vertical trans¬ 
mission is complete, then a unique endemic equilibrium exists and is not LAS when i?o < 1 and Rqu > 1. 
Backward bifurcation occurs when endemic equilibrium changes into disease free equilibrium if we decrease 
the initial number of infected individuals, or it reaches another LAS equilibrium if we increase the number 
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of initially infected individuals. When vertical transmission is incomplete but more than 50%, and < 1) 
two endemic equilibria coexist, but only one endemic equilibrium is LAS such that it becomes disease free 
equilibrium or endemic equilibrium by perturbation. Since Wolbachia-infected mosquitoes are less capable 
of transmitting dengue virus, complete Wolbachia-mfection is the ideal case for dengue control. 

When Rq > 1, Wolbachia can spread with a small number of initially infected mosquitoes, although 
very slowly. Population replacement may occur. When Rq <1, Wolbachia can spread if initial number of 
infected individuals exceeds a threshold, which depends on vertical transmission rate, ratio of egg laying 
rates of Wolbachia-mfected females to egg laying rates of Wolbachia-fiee female mosquitoes, and ratio 
of death rates of infected female mosquitoes to death rates of uninfected female mosquitoes. A smaller 
number of Wolbachia-iniected female mosquitoes is needed to be released for persistence of Wolbachia if a 
population suppression strategy is implemented before the release. 

The reproduction number is the product of the vertical transmission rate, ratio of the egg laying rates 
of Wolbachia-infected mosquitoes to egg laying rates of Wolbachia-free mosquitoes, and the ratio of death 
rates of Wolbachia-iree mosquitoes to death rates of Wolbachia infected mosquitoes. If the total number 
of eggs produced by Wolbachia-infected mosquitoes through vertical transmission is more than the total 
number of eggs produced by Wolbachia-fiee mosquitoes, then i2o > 1) such that complete infection or 
endemic equilibrium is LAS. 

The number of Wolbachia-iniected mosquitoes required for sustained Wolbachia infection depends on 
vertical transmission rate, the ratio of the number of eggs laid by Wolbachia-iniected mosquitoes to 
the number of eggs laid by Wolbachia-fvee mosquitoes, the ratio of the death rates of Wolbachia-free 
females to death rates of Wolbachia-infected females, and the ratio of the death rates of Wolbachia-fiee 
males to death rates of Wolbachia-infected males These parameters depend on pariculiar ITo/6ac/iia 

strain. If the life span of a mosquito is shorter, then the mosquitoes will lay fewer eggs. Once we know 
the specific Wolbachia strain that a specific mosquito species, such as Aedes aegypti or Aedes albopictus 
is carrying, we can estimate the number of infected individuals needed to be released for sustainable 
Wolbachia establishment using this model. 

We find that reducing both the aquatic stage and adult uninfected mosquitoes before releasing Wolbachia- 
infected female mosquitoes is the most effective strategy to reduce the number of Wolbachia-infected 
female mosquitoes needed for Wolbachia persistence (Table [6|). The second most effective strategy was to 
repeatedly kill the wild uninfected mosquitoes (to reduce both the adult and the aquatic stage mosquitoes) 
before releasing the infected mosquitoes. 

The model and analysis can help in understanding how Wolbachia can invade and persist in mosquito 
populations. In future research, we will couple our model to a dengue fever transmission model and analyze 
the impact that a bacteria infected mosquito population has on the spread of dengue. 
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Appendix 


The system of equations for mosquito population dynamics is: 

dF 

— = e-il^A - fifuF 

= (1 - 0)'ipA - HmuM 

The parameters for this model are described in Table [TJ 


(9a) 

(9b) 

(9c) 


Theorem 4.1 The zero equilibrium for mosquito population dynamics is LAS when Rqu < 1; while the 
steady state is LAS when Rqu > 1, where Rqu = • 


Proof. Jacobian matrix of system Q is: 


J = 


-Ha -fj- 
0lf 


1-^ 

Ra 

~Lfu 

0 


0 
0 

~Hmu_ 


The Jacobian matrix for no-infection equilibrium is: 

Jo = 


-Ha -fj 4> 0 

Ofj -Hfu 0 
(1 - 0)4! 0 -Hmu 


The characteristic polynomial of Jq is: 

(A -|- Hmu)[{^ + HaP V')(A + Hfu) ~ 4>^4^] = 0 

If Rqu < 1, then all eigenvalues are negative, the zero equilibrium is LAS. 
The Jacobian matrix at steady state is: 


Jss - 


-Ha -4^- 4 >t^ 4)[i- 

04; -Hfu 

{1-0)4; 0 


0 
0 


Characteristic polynomial of Jss is: 

{)^ + Hmu) X + Ha + 4^ + (l- 1 {X +Hfu) - (t>04^]=^ 

[ Hfu V yj 

If Rqu > 1, then all eigenvalues of Jss are negative, and the non-zero steady state is LAS. 


( 10 ) 


( 11 ) 
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Theorem 4.2 The zero equilibrium: (0,0,0) is globally asymptotieally stable (GAS) when Rqu < 1- The 
steady state: {Ka{l - 7 ^), 7^(1 - 7 ^), “ 7 ^)) GAS when Rou > 1- 


Proof. We follow the approach in 3^. First we consider the subsystem: 


dt 

dt 


= 0F 1 - 


A 


Ka 

= Of;A - fifuF 


- piaA - 'IpA 


( 12 a) 

( 12 b) 


where / = {A, F) is on an open set D C 

Let (5i = ^2 = Ij 1^ 1 ^ = Of: > 0, = </> > 0, and (^ 2 ^ =///„> 0. According to the definition of 

tridiagonal feedback 3^ , system [12] is a monotone tridiagonal feedback system with Poincare-Bendixson 
property [37 1. 


Recall Theorem 2 in 


38l |. If the systems of ODEs dx/dt = f{x),x C D satisfies: 


(1) The system exists on a compact absorbing set K d D. 

(2) A unique equilibrium point E exists and is LAS. 

(3) The system has Poincare-Bendixson property. 

(4) Each periodic orbit of the system is asymptotically stable. 

Then E is globally asymptotically stable in D. 

To prove that each periodic orbit Q = p{t) : 0 <t < w of system (fT^ is asymptotically stable, we follow 
39 1 and Theorem 3 in [^. We need to prove that the linear system = j'p{p{t))z{t) is asymptotically 


\ 2 ] 

stable, where Jp is the second additive compound matrix of the Jacobian matrix Jp associated with 
system m- Eor system m, = — (^4> + P'a + '^F + . We build the following linear system with 

one equation and the right hand side is the compound matrix of the Jacobian matrix Jp. 

— - - {'Ip + tia + —F + Hfu ] X. 


Kr, 


Let Lyapunov function V(X, A,Y) = |X|. 

The right derivative of P along the solution paths (X) and {A, F) is: Dp{V{t) = —{^p+|aa+^F+f^fu)\X\, 
which implies that V{t) 0, and A(f) ^ 0 as t —)• 00 . Therefore, the linear system [12] is asymptotically 
stable, and the solution (A, F) is asymptotically orbitally stable with asymptotic phase. 

By the same argument in 3^, the system [12] is uniformly persistent in C The zero equilibrium 
(0,0,0) is isolated and the largest compact invariant outside D is (A*.F*), which is absorbing and the 
system is uniformly persistent [^. The conditions for Theorem 2 in |38[ | are all satisfied. Therefore, 
(0,0) is GAS whenever Rqm < 1; and (A*, F*) exists and is GAS when Rq^ > 1- 

Following ^1 and Theorem 4 in 36], (0,0,0) is GAS whenever Rqu < 1 and {A*,F*,M*) is GAS 
whenever Rqu > 1 - * 


Theorem 4.3 



/ 



0 < A < \ 

II 

_/s_ 

< 

F 

\m) 

G 

0 < F < , 

0 < M < 1 

fAmu ' 


( 13 ) 
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is an invariant region under the flow induced by m- 


Proof. The proof directly follows the proofs for Lemma 4.2 and 4.3 in 421. 
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